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In this paper we analyze the biasing effect of point sources, either thermal Sunyaev-Zeldovich clus- 
ters or standard radio sources, on the estimated strength of the non-Gaussianity in the Cosmic Mi- 
crowave Background (CMB). We show that the biggest contribution comes from the cross-correlation 
... of the CMB with the matter density rather than from the poisson term which is conventionally as- 

00 ■ sumed in these calculations. For the three year WMAP data, we estimate that point sources could 
' produce a non-Gaussian signature equivalent to a bias in /jvl of 0.35,0.24, —0.097, —0.13 in the 
' Ka, Q, V and W bands respectively. The level of bias we find is largely insufficient to explain 

the very high /atl values recently detected by Yadav and Wandelt. For Planck, we estimate the 
5-H ' point source bispectra to contaminate the /jvl estimator with a bias of 1.3,0.34, —0.25, —0.48 at 

^2 \ 30, 44, 70, 100 GHz respectively. These results depend on the assumed redshift distribution of the 

. point sources. However, given the projected Planck sensitivity of A/jvl ~ 5 (95 % C.L.), a good 

' estimate of point sources' properties including their number density and redshift distribution is 

essential before deriving strong conclusions on primordial non-Gaussianity. 

<N : 

. PACS numbers: 98.70.Dk, 98.80.-k 

6 ■ 

^ ■ I. INTRODUCTION 

C/3 ! 

1 Recent claims by Yadav and Wandelt [Ij of the detection of strong primordial non-Gaussianity in the three- year 
Wilkinson Microwave Anisotropy Probe (WMAP) data Q have the potential to revolutionize our understanding of 
the early universe. These results were also found in the WMAP five-year analysis, although with less statistical 
significance [29]. The strength of the non-Gaussianity detected in their analysis is more than two orders of magnitude 
larger than the non-Gaussianity expected in the simplest model of single field, slow roll model of inflation [3[. This 
detection, if it stands up to scrutiny, will be the first definitive indication that the simplest model of inflation cannot 
adequately explain all of the current cosmological observations and must be modified in some way. 

In Yadav and Wandelt's analysis the detection of the non-Gaussian signal is due to the simultaneous reduction 
^— ^ in the estimator's error bars and the shift in the central value as smaller angular scale information was included in 
the analysis. One obvious concern is that the estimator is contaminated by foreground emission, in particular radio 
, points sources and the thermal Sunyaev-Zeldovich effect which become increasingly more important on small angular 
' scales. 

In this paper we will analyze the influence of point sources, for the experiment resolution of WMAP both the radio 
' sources and SZ clusters effectively act as point sources, on the standard non-Gaussianity estimator. While it has been 
Jh , claimed that the non-Gaussianity caused by Poisson fluctuations in the number density of radio point sources can be 
. . I safely separated from the primordial non-Gaussian signal [3], we will demonstrate the other forms of point source non- 
Gaussianity cannot be safely ignored. In addition to the standard forms of non-Gaussianity produced by point sources, 
we will show that cross-correlation between the point source power spectrum and the CMB temperature anisotropies 
or instrument noise will produce non-Gaussianity of a very similar form as the local model. The local model is the 
form typically assumed in analyses of primordial non-Gaussianity, so these new non-Gaussian contributions will bias 
the estimator in a fashion that cannot easily be corrected. 

This paper is organized as follows. In Wl\ we demonstrate that point sources can produce a bispectrum that has 
the same form as the local model. In ^jll]we derive the bias induced by the various point source bispectra. In WII we 
conclude. We use the WMAP three-year cosmological model 01 for numerical calculations. 
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II. ESTIMATOR BIAS 



While non-Gaussianity generically implies that any higher-order connected correlation function is non-zero, it is 
typical to focus on the three-point correlation function, or equivalently the bispectrum, because it has the simplest form 
of all non-Gaussian correlation functions and for weak non-Gaussianity it contains the nearly all of the information 
The three-point correlation function can be factored into a component fixed by rotational invariance, which is assumed 
a prioiri, and a piece determined by the underlying mechanism that produced the non-Gaussianity [4|. Rotational 
invariance forces the three-point correlation function to be proportional to the Gaunt integral. 
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and the CMB three-point correlation function can be written as 



(afimiaf2m2««3m3) " ^mimsmg 1 ,^2 A 



(1) 



(2) 



where the reduced bispectrum, 2,^31 contains information about the form of non-Gaussianity. 

The estimators used in CMB non-Gaussianity analyses are optimized for the detection of a signal with a very 
particular form, namely the local model [H]- The local model assumes that the initial curvature perturbations can be 
written as 



^{x)=C^g{x)+fNL[<l^l{x)~{4>l{x))]. 



(3) 



where (pg is Gaussian. The non-linear terms in this model lead to the following bispectrum for the initial curvature 
perturbations 



B{ki,k2,k3) = 2fNL[P{ki)P{k2) + eye]. 



(4) 



where P{k) is the power spectrum. The ordinary linear radiative transfer function are subsequently used to calculate 
the CMB temperature anisotropics from these initial curvature perturbations. The statistical properties of the CMB 
temperature anisotropics will mirror the statistical properties of underlying curvature perturbations since we are 
consider linear radiative transfer. The levels of non-Gaussianity claimed by Yadav & Wandelt are significantly larger 
than the expected non-Gaussianity produced by non-linear radiative transfer. 

However the estimators are sensitive to any bispectrum form that might be present in the data, regardless of its 
origin. In this section we will determine the induced bias produced by the various cross-correlation bispectra described 
in the next section. 

The non-Gaussianity estimator can be expressed as @ 



fNL-J 



where the normalization is 
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(6) 



here Cj = Ci + is the sum of the CMB signal and noise. The CMB experimental noise parameters are described 
in Table U 

We are ignoring the additional linear term in the estimator because the contribution of radio point source will not 
bias this piece of the estimator if the signal and noise covariance matrices well represent the real data. 

The weight functions used in the estimator are optimized for a bispectrum produced by the local model. The 
estimator bias will be determined by substituting the various forms of the point sources bispectra into the estimator 
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Here ^t^^i^^i^ is the one of the reduced bispectrum produced by point sources. The CMB bispectrum produced by the 
local model is generally negative. The collapsed triangle modes, which have the highest signal-to-noise, are always 
negative. So a positive point source bispectrum will cause a negative estimator bias. We will now discuss the possible 
bispectrum terms. 
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III. POINT SOURCE BISPECTRA 

The observed signal is the sum of the primordial and secondary temperature anisotropies, foreground emission and 
instrument noise. While the secondary anisotropies and extra-galactic foregrounds, which maybe quite non-Gaussian, 
are important on small angular scales, the signal on the large angular scales relevant for WMAP is dominated by 
primary anisotropies. Thus it is assumed that any measured non-Gaussianity by WMAP is primordial in nature. We 
will argue that cross-correlations between some of these signals may induce bispectra on large angular scales. 

Moreover these bispectra may have similar forms as the local model if the signal power spectrum becomes spatial 
inhomogeneous in manner that then correlates with a second component in the observed signal. For example, the 
matter overdensity will bias the local radio point source power spectrum and it will be correlated with the CMB 
temperature anisotropies produced via the ISW effect. This will produce a bispectrum similar in form to the local 
model. 

This is not a coincidence as the non-Gaussianity in the local model is produced by the modulation of the small 
scale inflaton power spectrum by the large scale inflaton fluctuations that have already left the horizon and frozen 
out. The parameter /^vl is a measure of the non-linear coupling between these different scales. Likewise the large 
scale matter overdensity modulates the small scale Poisson fluctuation power spectrum by altering the local number 
density of radio point sources. In an analogous fashion the bias describes the coupling between the different scales. 
As discussed in the previous section any bispectrum present in the data can bias the estimator. 

IV. ANISOTROPY MECHANISMS 

Now we will discuss the various physical effects considered in this paper - radio point source emission, the thermal 
Sunyaev-Zeldovich effect and the integrated Sachs- Wolfe effect. 

1. Radio Point Sources 

The temperature anisotropy induced by unresolved radio point sources can be expressed as an integral over their 
flux distribution function 

AT 1 /■■^^"'^ dN 

— {n,y) = - dSS—{S,v-h). (8) 

T Cy Jq dS 

The conversion between the temperature and intensity fluctuations is 

where Tcmb — 2.728 K and B^{T) is a blackbody frequency distribution. 

Note that we have allowed both the upper flux limit for unresolved radio point sources and their number density to 
be spatially inhomogeneous. As we will describe below, these spatial inhomogcneities will correlate with other signals 
present in the data to produce bispectra in the observed CMB data. 

The actual values for the bispectra calculated in this paper depend on the radio point sources properties (both flux 
and redshift distributions) and on the flux cut at a given frequency for a specific experiment. In order to evaluate the 
residual point source contribution to the total estimated bispectrum in the WMAP-Syrs data, we must consider the 
technique used by WMAP to identify and substract radio point sources and the point source mask applied to the 
data. The WMAP source selection criterion does not correspond to a single flux threshold at a given frequency, rather 
it requires that a candidate source should be seen with a minimal statistical signiflcance in all channels Q . As Yadav 
& Wandelt's results [Tl, which are the motivation for this paper, were derived using the V and W frequencies bands, 
the radio point source populations at these frequencies is the most relevant. Unfortunately, as most radio sources are 
stronger at low frequencies, they tend to be detected with higher significance in the K - Q bands than in the V - W 
ones. 

In addition from a blind search, at 20-30 GHz it is possible to use lower-frequency catalogs of point sources as 
tracers for detection ^. As a result, the source population at such frequencies is much better characterized than the 
one at higher frequencies. By considering fiux number counts at low frequencies it is possible to give an estimate of 
the fiux above which the detected point sources create a complete catalog. This flux threshold is estimated to be 
above 1.1 Jy at 23 GHz (K Band) [9|, HO)], while in the W band (94 GHz) the number counts are not sufficiently well 
determined to allow such estimate. An alternative blind search technique applied to the WMAP V and W band data 
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Planck - 30 


0.33 


33 
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Planck - 44 


0.36 


23 


2.4 


Planck - 70 


0.34 


14 


3.6 


Planck - 100 


0.13 


11 


1.6 



TABLE I: Radio point source flux threshold, FWHM, and instrument pixel noise (in 10 ) for the relevant WMAP and Planck 
frequency bands. For WMAP the various band names are also listed. 

increases the number of sources found in the V band by 50% compared to the WMAP team results [TT| . These results 
inchide some sources that were not contained in the WMAP point source mask. The new sources found, however, do 
not seem to represent a different population than the ones previously detected. As the WMAP 5 year data is now 
available, more work on point source characterization at the frequencies where the CMB science is derived should be 
possible. 

The actual residual contribution of point sources depends upon the mask applied to the observed map. In the case 
of WMAP, this mask considers the WMAP detected sources (about 300) as well as some bright sources from other 
low-frequencies catalogs, for a total of approximately seven hundred sources masked. The actual selection function 
that this procedure imposes at the V and W bands is poorly understood, and indeed some sources detected by [ll| 
were not masked. However, while the detection threshold implied by this whole procedure is poorly defined, there is 
good agreement in the residual point sources power spectrum contribution as derived by different authors 0, [ill , [l^ . 
This can be approximately translated in a flux threshold of 0.6 Jy in the Q band 12] and 0.75 Jy in the V band 
For illustrative purposes, we will adopt here an approximate estimate for the flux cut-off of 0.7 Jy in all WMAP 
bands. In f|V]we discuss the dependence of our results on this choice. 

The Planck satellite will have better resolution and sensitivity, resulting in a lower detection threshold for point 
sources. In the following, we make predictions of the bispectrum expected in the final Planck maps, considering 
the detection thresholds for a 95% complete sample derived by jl^l using realistic Planck sky simulations. The flux 
cut-offs for both WMAP and Planck, as well as, the adopted instrument noise parameters are given in Table HI 

Finally, in order to compute the bispectra implied by point sources below a given flux, we adopt the source counts 
predictions of Toffolatti et al [l3| rescaled by a factor 0.8, as suggested by the matching of these predictions with 
the actual number counts at fluxes above 1.1 at 41 GHz obtained by fld\. We will use the radio point source bias 
f^PS 2^ J [isl. [igI, [17| as inferred for low frequency radio point sources. This value is uncertain and almost definite 
varies for the various population type that constitute the high frequency sample. 

It is difficult to constrain the redshift distribution due to the lack of optical studies of 60 — 96GHz source population, 
work at 23GHz has been conducted by Gonzalez et al. for fluxes above 1 Jy [lo|. At these fluxes and frequencies 
the number counts are dominated by QSOs for relatively high redshifts. This population is well approximated by the 
following analytical formula: 

n^^(z) oc 0.75 X e-(--o)V(2-)^ ^ (10) 

with zq = 0.95 and a is 0.4 (0.9) for z < 1 (z > 1). This redshift distribution is in agreement with recent studies of 
radio sources populations at 90 GHz with ATCA In addition, there is a small (10-15 %) of the total population 
that consists of radio loud galaxies which peaks at much lower redshifts (z < 0.1). Given the small number of sources, 
it is difficult to derive an appropriate fitting formula for this other population. In the following, we will adopt the 
following distribution of low redshift sources 

n^^(z) cx 0.25 X 10-^^ (11) 

that provides a better fit to counts found by Gonzalez et al. [l3|- We will take the sum of EqslTUl&IIIlwith a common 
proportionality factor determined by requiring n^^{z) to integrate to unity over the range z = to z = 3.1. The 
relative amplitude of the QSO and radio galaxy contribution to the sources is derived from the optical indentifications 
of Gonzalez et al [l3|. This fit will be called model 1. Gonzalez et al. [l^|_provide a fit to a theoretical model of the 
luminosity functions of the various source populations as derived from l20l|. To understand how the uncertainty in 
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the source redshift distribution affects our results, we also do calculations with the model of Gonzalez [lOj; this will 
be called model 2. 

WMAP only resolves the high- flux sources which are typically dominated by AGN. Most likely lower flux sources 
consist of a different population (e.g. [T^, '20'|) and therefore have a different redshift distribution with possibly more 
weight either at lower or higher redshifts. An analogous redshift analysis on higher frequencies catalog is strongly 
needed, but is not available at this time. Future investigations of radio point sources catalogs plus WMAP and Planck 
results are higher frqeuencies will help clarify this issue. For the aims of this paper we take Eq. [TOlto be the redshift 
distribution at all frequencies and fluxes and keep in mind the potential uncertainty that this assumption introduces. 



2. Thermal Sunyaev-Zeldovich Effect 

The hot plasma in the intra-cluster medium will produce temperature anisotropics via Thomson scattering of the 
incident CMB photons; this is the well-known thermal Sunyaev-Zeldovich (SZ) effect (see Carlstrom et al. [2l| for a 
review). Following the model of Komatsu & Seljak [13, [1^, the temperature anisotropics produced by the SZ effect 
can be expressed as an integral over the cluster mass distribution function 

— {n,iy)^g,Jdz—j^ dMy{M, z)^[M, z;nx{z)], (12) 

here the frequency dependence of the thermal SZ effect is given by 

e"^ 4- 1 

gu = x^—-^, (13) 
— 1 

where x = hv /ksTcus,- The Compton y-parameter is related to the line-of-sight integral of the the cluster's thermal 
pressure, dn/dM is the cluster mass function and the volume clement is 

-d^ - im^^'^^ 

where x(^) is comoving distance to redshift z. The details of the implementation of this model are extensively discussed 
in (2^ . Ignoring clustering terms in the power spectrum, we find the thermal SZ point source power spectrum 

dz—J^ dMyj{M,z)^{M,z). (15) 

In a similar fashion to the radio point sources, the number density of massive clusters will be changed by both the 
biasing effect of the large scale matter overdensity and gravitational lensing magnification. In order to determine how 
these processes will affect the SZ power spectrum we need to know the redshift distribution of power in the SZ effect. 
The weighted redshift distribution of power produced by the SZ effect can be expressed as the integral over halo mass 
of the Sheth-Tormen halo mass function 

n'^^(z) oc — / dMM^" — {M,z). (16) 

dZ J Myriia 

The cluster y-parameter - mass scaling relationship slope is taken to be a = 1.6 jHf. The limits of integration are 
taken to be Mmin — 10^'* M0 and Mmax = 5 x IQ^^Mq. The redshift distribution is normalized so it will integrate to 
unity. We also need the bias weighted redshift distribution of power produced by the SZ effect can be expressed as 
the integral over the Sheth-Tormen halo mass function 

{bnf^{z)oc^ [ "^"^^ dMM^''-$^{M,z) b{M,z), (17) 

dz JMmin dM 

where b{M, z) is the standard bias of the Sheth-Tormen mass function. 



3. Integrated Sachs-Wolfe Effect 

Most of the power in the CMB temperature anisotropics is produced at high redshift during recombination. There 
will be very little cross-correlation between these temperature anisotropics and the low redshift matter overdensity re- 
sponsible for altering the small scale point source power spectrum. However additional CMB temperature anisotropics 
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FIG. 1; Redshift weight functions for the ISW effect (green, dot-dashed); radio point sources - model 1 (red, dashed) and 
model 2 (blue, long-dashed); and thermal SZ effect (solid, black). 

can be generated at low redshift via the Integrated Sachs- Wolfe (ISW) effect if the gravitational potential fluctuations 
are evolving in time. In a fully matter dominated regime the gravitational potential fluctuations are static, however 
as the universe becomes dark energy dominated the ISW effect can occur. 
The temperature anisotropy produced by the ISW effect can be expressed as 

-jrin) ^ -2j dz — , (18) 

= iHl^lM j dz[{l + z)D{z)]' j l^^e^^-^'^^, (19) 

where D{z) is the linear theory growth function and the prime denotes differentiation with respect to z. 

In Fig. [l]we show the redshift weight functions for the ISW effect [(1 -I- z)D{z)]' (green, dot-dashed); the radio 
point sources b^^n^^ {z)D{z) - model 1 (red, dashed) and model 2 (blue, long-dashed); and the thermal SZ effect 
(bn)^^ {z)D{z) (solid, black). The overlap of these weight fuctions will determine the amplitude of the cross-correlation 
spectra as described in i^IVBI fc flVCI and shown in Figs. [1]&|31 

A. Point Source Bispectrum 

The simplest bispectrum form produced in the CMB is due to radio point source Poisson fluctuations. The reduced 
bispectrum is independent of scale and can be written as 

KM,t,=cZ^ dSS' — {S,u). (20) 

This bispectrum component has been detected in the WMAP data 26'] . Since its functional form is not similar to the 
local model it will not significantly bias the estimator and we will ignore it as did Yadav & Wandelt. The five-year 
WMAP analysis [2^ includes estimates (and corrections) for the estimator bias produced by this bispectrum form. 
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B. Number Density Modulation 

The Poisson fluctuation power spectrum in a certain region of the sky can be written as integral over the distribution 
of sources in that region. If the anisotropic component of the power spectrum correlates with any other signal present 
in the data, non-Gaussian correlation functions will be generated. In this subsection we will focus on correlation of 
the large scale matter overdensity, which bias the number density of point sources, with the ISW effect. 

The power spectrum produced by radio point sources in direction h is 



where the anisotropic point source distribution can be expressed in terms of the matter overdensity as 



dN dN 

— iS,.;n)^—iS,.) 



dzn^^{z)6{n, z) 



(21) 



(22) 



the mean point source distribution was described in ijllll The large scale matter overdensity, which biases the local 
number density of point sources as shown in Eq. [211 will be correlated with the large scale CMB temperature 
anisotropics produced by the ISW effect and the following reduced bispectrum will be induced 



(23) 



Here the isotropic source distribution dN/dS leads to an isotropic power spectrum C^^ . The matter-ISW cross 
correlation spectrum can be expressed as 



X 



PS 



PS 



H{z)dz 



P 



X{z) 



D{z)n''^{z)[{\ + z)D{z)]' 



(24) 



where P{k) is the matter power spectrum and h^^ ~ 1.7 is the radio point source bias. The cross-correlation spectrum 
is positive because matter overdensities correspond to potential wells. At low redshift as the amplitude of the potential 
wells decay the CMB photons experience a net blue-shift, thus the cross-correlation is positive. In Eq. [21] we have 
employed Limber's approximation to simplify the cross-correlation spectrum. The function n^^{z) is the redshift 
probability distribution function of the radio point sources defined in Eq. [TOl The large scale matter overdensity 
affects the thermal SZ power spectrum in an analogous fashion. 

In Fig. [21 we show the matter-ISW cross correlation spectrum for the thermal SZ effect (solid, black) and the radio 
point sources - model 1 (red, dashed), model 2 (blue, long-dashed). The thermal SZ effect has the largest cross 
correlation because the clusters tend to be located at lower redshift than the radio point sources. The ISW effect 
primarily occurs at low redshift once the universe is strongly dark energy dominated, so it most strongly overlaps 
with the SZ effect. 

In Fig. [3] we show the number density modulation bispectrum for the thermal SZ effect (solid, black) and the radio 
point sources - model 1 (red, dashed), model 2 (blue, long-dashed). Also shown for reference is the equilateral shape 
(£i = £2 = ^3 = ^) of the primordial local model bispectrum {/nl = 1) (green, dot-dashed). The equilateral shape of 
the local model bispectrum changes signs, the zero-crossing are obvious from the plot. This is a consequence of the 
radiative transfer functions producing both hot and cold regions on the sky. The collapsed shape always has the same 
sign, opposite the sign of /nl- For the local model the collapsed bispectra have the highest signal-to-noise. 



C. Magnification Modulation 



The distribution of matter along the line-of-sight between the observer and the sources will gravitationally lens 
these point sources. The priniciple effect of gravitational lensing will be to change the source density by magnifying 
and de-magnifying certain regions of the sky. The magnification in a given direction can be written, to first order, as 



/x(n) ~ 1 + 2k,{ji). 

The convergence field distorting a background source at z is related to matter overdensity as 



K(n, z) 



,n,,Hl /|f_(i + .0^[x(.)-x(.')]^(n,.'), 



(25) 



(26) 
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FIG. 2: Density-ISW cross-correlation spectrum for thermal SZ (solid, black); radio point sources - model 1 (red, dashed) and 
model 2 (blue, long-dashed). 




1 



FIG. 3: Number density modulation bispectrum for thermal SZ (solid, black); radio point sources - model 1 (red, dashed) and 
model 2 (blue, long-dashed). The equilateral shape {£i = £2 = is = i) of the local model (/jvl = 1) (green, dot-dashed) is also 
shown for reference. 
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and the average magnification of SZ point sources in a given direction is then given by 

.(n)./d.n-(.Mn,.). 



The observed SZ cluster number density in a given direction will be related to the intrinsic number density as 

dn 



iM,n), 



dn 



^ [l-2K(n)] — (M,n). 

The convergence field is correlated with the CMB temperature anisotropics produced via the ISW effect, 



SZ 



2^2 



dzn^^{z) / dziP 

"'0 



Xizi) 



{1 + zi)D{zi)[{l + zi)D{zi)] 



/X(£i) 
Xiz) 



[x{z) - x{zi)]- 



(27) 

(28) 
(29) 

(30) 



Again we have evaluated the cross-correlation according to Limber's approximation and the convergence-ISW cross- 
correlation spectrum is positive. This magnification effect will result in the following reduced bispectrum 



-ySZ 



SZ (riSZ 



-ySZ\ i,,rSZ ir<SZ 



',SZ\ 



(31) 



The bispectrum is negative because large scale CMB hot spots (positive ISW effect) correlate with positive magnifi- 
cation which always reduces the amplitude of point source Poisson fluctuations. 

The radio point source selection function is expressed in terms of the observed flux which can be affected by 
gravitational lensing due matter along the line-of-sight. The gravitational lensing magnification will modulate the 
flux cut-off and be correlated with the ISW temperature anisotropics. The point source dilution effect discussed above 
will also occur. The radio point source Poisson fluctuation power spectrum can be expressed as 



dSS^ 



1 dN 



Linearizing in the convergence field, we find that the Poisson fluctuation power spectrum becomes anisotropic 



This will result in the following reduced bispectrum 



-4(Mr 



eye.) 



(32) 



(33) 



(34) 



where Mf^ is the radio point source version of Eq. [30l Note that gravitational lensing affects the radio point source 
power spectrum by changing both the upper flux cut-off and the source counts, whereas the SZ power spectrum is only 
altered by changes in the local cluster counts. If the SZ clusters are detected with high signal-to-noise and removed 
from the CMB maps then flux cut-off modulation effect will also produce an additional bispectrum term. 

In Fig. Uwe show the convergence-ISW cross correlation spectrum for the thermal SZ effect (solid, black) and the 
radio point sources - model 1 (red, dashed), model 2 (blue, long-dashed). The model 2 of the radio point sources 
has the largest cross correlation because it predicts that the point sources tend to be located at higher redshift. As 
opposed to the matter-ISW cross-correlation, which requires that the sources lie in the same redshift range over which 
the ISW effect occurs, in this case the point source do not have to be at the same redshift at which the ISW occurs 
in order for the contribution to be relevant. In fact the point sources experience greater magnification if they are at 
substantially higher redshift than the matter distribution that is simultaneously magnifying them and is correlated 
with the ISW effect. 

In Fig. [5] we show the magnification modulation bispectrum for the thermal SZ effect (solid, black); radio point 
sources - model 1 (red, dashed) and model 2 (blue, long-dashed). The local model bispectrum {/nl = 1) (green, 
dot-dashed) is also shown for reference. 
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FIG. 4: Convergence-ISW cross-correlation spectrum for thermal SZ (solid, black); radio point sources - model 1 (red, dashed) 
and model 2 (blue, long dashed) 




FIG. 5: Magnification modulation bispectrum for thermal SZ effect (solid, black); radio point sources - 
and model 2 (blue, long-dashed). The local model bispectrum {/nl = 1) - equilateral shape {£i = 
dot-dashed) is also shown for reference. 



model 1 (red, dashed) 
^2 = ^3 = i) (green. 
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D. Selection Modulation 



The amplitude of the Poisson fluctuation power spectrum depends on the number density of radio point sources 
below some flux Umit determined by the radio point source removal technique. If the selection criterion used produces 
an anisotropic flux limit that is correlated with either the large scale temperature anisotropy or the instrument noise, 
then a bispectrum will be produced. If the radio point sources are identified via external catalogs, such as NVSS or 
PMN, then there will be no cross-correlation and therefore no bispectrum. And if multi-wavelength data is differenced 
to isolate the power-law frequency dependence of the radio point source, then the modulated selection function will 
produce correlations between the data in different frequency bands. 

The simplest method is to removal pixels above some multiple (7) of the pixel variance 



S{n) = Ci, 



/AT 
V T 



(n) + N{n) 



(35) 



The temperature anisotropy and instrument noise have vanishing expectation values, so the mean flux cut-off is 

So = c„j(Tp, (36) 

where the pixel variance is deflned as 



= [Any y: 



An 



C. 



(37) 



here Cj = + is the sum of the CMB signal and noise and AO, is the pixel size. There will be fluctuations 
about the expected flux cut-off that are correlated with either the temperature anisotropy or the instrument noise. 
Linearizing in these fluctuations, we find the following reduced bispectrum 



-^{Cl+cyc.)-sf-^{S,u). 



(38) 



The selection criterion used by WMAP is based on an algorithm developed by Tegmark & de Oliveira-Costa 
[27I I . The map is filtered in order to reduce the importance of the long- wavelength CMB modes. The filtered total 
temperature fluctuation in some pixel is 



y 



(n) = ^ yfpYim{fi){atr^ 



)An + J2 



(2^+1) S{n) 



(39) 



The algorithm removes any pixel that has a value greater than some multiple (7) of the flltered map pixel variance 

(2^+1) 1 



af = {An)-^Y.- 



47r Cj 



(40) 



The threshold 7 is chosen according to some compromise between false positive and negatives. Since the thresholding 
is applied to the total signal in a pixel, the corresponding radio point source flux cut-off in a pixel is 



S{h) 



F 



l<Jp 



a£m + nim) 



(41) 



where the normalization is 



(21+1) 1 



, 4. CJ- 

The temperature anisotropy and instrument noise have vanishing expectation values, so the mean flux cut-off is 

So = y70-p- 



(42) 



(43) 



There will be fluctuations about this expected value and these fluctuations will be correlated with either the temper- 
ature anisotropy or the instrument noise. Linearizing in these fluctuations, we flnd the following reduced bispectrum 



3An -^dN - 
c^b dS 



(44) 
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In regions with large instrument noise or CMB temperature anisotropy, the radio point source flux cut-ofF will be 
lowered. This reduces the total number density of radio point sources in that region and therefore the Poisson 
fluctuation power spectrum. This effect explains the negative sign in the reduced bispectrum, Eq. O 

Since the filtering applied to the maps produces a bispectrum independent of scale, similar to the point source 
bispectrum described in i]IV A[ we will ignore it. In reality the actual selection function is more complicate than this 
simple filter technique predicts and it might produce a bispectrum of a much different form. Numerical simulations 
incorporating the exact selection procedure will need to be done in order to fully determine its effect on the estimator. 

V. NUMERICAL RESULTS 

In this section we will present numerical results for the contamination of the standard non-Gaussianity estimator 
by the various bispectra discussed in this paper. This is done for both the WMAP and Planck instrument noise levels, 
frequency bands and flux cut-offs given in Table [D 

In Fig. [5]the f^^ estimator bias as a function of imax is shown for the different bispectra - radio point source number 
density modulation (solid, black); SZ number density modulation (dotted, red); radio point source gravitational lensing 
magnification modulation (dashed, blue); and SZ gravitational lensing magnification modulation (long-dashed, green) 
with WMAP instrument noise. The bias plots are shown for the four relevant WMAP frequency bands - upper left 
Ka - 33 GHz; upper right Q - 40 GHz; lower left V - 61 GHz; lower right W - 94 GHz. The magnification modulation 
effect produces a positive bias since its bispectrum is negative, while the density modulation effect produces a negative 
bias. 

Yadav & Wandelt claim a central value of Jnl — 86.8 with a standard deviation of cr = 30.0 for a 2.9a detection of 
non-Gaussianity [l[. In the five-year WMAP data the central value is found to he Jnl = 67 with a standard deviation 
of cr = 31 [1^. At imax = 750 the total estimator bias A/atl = 0.35 in the Ka band, AJnl = 0.24 in the Q band, 
A/nl = —0.097 in the V band and AJnl = —0.13 in W band. Since the density modulation and the magnification 
modulation bispectra have different signs they partially cancel and reduce the overall effect. At low frequency the 
radio point source magnification modulation bispectrum is the most important so the bias is positive. At higher 
frequencies the SZ density modulation bispectrum dominates which makes the bias negative. These numbers should 
be compared to the estimator bias produced by the radio point source Poisson fluctuation bispectrum. Komatsu et 
al [i^ have estimated this bias to be AJnl ~ — (3 — 5) at imax = 700. 

In Fig. [7]the fj^^ estimator bias as a function of imax is shown for the different bispectra - radio point source number 
density modulation (solid, black); SZ number density modulation (dotted, red); radio point source gravitational lensing 
magnification modulation (dashed, blue); and SZ gravitational lensing magnification modulation (long-dashed, green) 
for Planck at 30 GHz (upper left); 44 GHz (upper left); 70 GHz (lower left) and 100 GHz (lower right). 

Estimates of the sensitivity on Jnl achievable with Planck suggest that A/jvl — 10 at 95 % C.L. using just 
temperature information and A/jvl ~ 5 at 95 % C.L also including polarization. Summing the various biases we find 
A/atl = 1.3 at = 30 GHz, AJnl = 0.34 at ;/ = 44 GHz, AJnl = -0.25 at i/ = 70 GHz and AJnl = -0.48 at 
V = 100 GHz. These results imply that a good knowledge of point source properties is important if Planck will be 
able to achieve its full potential in constraining primordial non-Gaussianity. 

There are uncertainties in the models we have used to describe the radio point sources. These model uncertainities 
will directly lead to uncertainities in the above predictions. As can be directly seen in Figs. [2]-[5l the differences 
between the radio point source redshift distributions, model 1 and model 2, are not significant. This is not surprising as 
the ISW kernel, with which the radio point source redshift distributions are being cross-correlated, is quite broad. The 
redshift distributions are supposed to trace both the high-fiux source populations at low frequencies. However, they 
may not be representative of the source populations at all frequencies considered here for lower flux cut thresholds. If 
a lower flux cut implies a higher population of low-redshift objects, the amplitude of the ISW-density cross-correlation 
spectrum would be increased. The decrease in the flux cut-off will decrease the Poisson fluctuation power spectrum, 
so the change in the bias of the /nl estimator is not clear. We also note that reducing the flux cut-off from the 
WMAP to the Planck level does not reduce the bias implied on /atl as the Planck noise levels and beam sizes are 
also smaller. 



VI. CONCLUSION 

In this paper we analyzed the effect of point sources, both due to radio emission and the cluster SZ effect, on 
the estimation of primordial non-Gaussianity in the cosmic microwave background. The standard non-Gaussianity 
estimator is sensitive to any bispectrum present in the data. In addition to the standard Poisson fluctuation bispec- 
trum, we found that cross-correlations between the radio point source and SZ power spectra and either the CMB 
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FIG. 6; The WMAP estimator bias terms A/^^ as a function of imax for radio point source density modulation (solid, black); 
SZ number density modulation (dotted, red); radio point source gravitational lensing magnification modulation (dashed, blue); 
and SZ gravitational lensing magnification modulation (long-dashed, green) in Ka - 33 GHz (upper left); Q - 40 GHz (upper 
right); V - 61 GHz; and W - 94 GHz. The density modulation terms produce a negative bias, while the magnification bias 
produce a positive bias. 
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FIG. 7: The Planck estimator bias terms A/j^^ for 30 GHz (upper left); 44 GHz (upper left); 70 GHz (lower left) and 100 GHz 
(lower right). The curves are the same as Fig. |6] The density modulation terms produce a negative bias, while the magnification 
bias produce a positive bias. 
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temperature anisotropies or instrument noise can produce bispectra. These bispectra have forms somewhat similar 
to the local model, which is the standard bispectrum form used to search for primordial non-Gaussianity in CMB 
data. These similarities are not accidental, but occur because the same basic principle generates the non-Gaussianity 
in these cases. Due to this similarity it will be much more difficult to distinguish this non-Gaussianity from the 
primordial signal than other secondar y b ispectra with shapes which can be quite different. 



A related paper by Serra & Cooray [28| has examined different secondary bispectra, but has reached similar conclu- 
sions. They examined the bispectra produced by the cross-correlation of the thermal SZ effect and the gravitational 
lensing of the primary CMB anisotropies. They concluded that the effects are too small to account for the infered 
values of Jnl from the WMAP data and will start to become important for the Planck dataset. 

The estimator bias that we have calculated is quite small and is not able to explain the results found by Yadav & 
Wandelt, although the bispectra considered in our paper will start to become important for the Planck non-Gaussianity 
analysis. Due to the extremely important nature of the detection of primordial non-Gaussianity in the WMAP data, 
all possible alternative explanations should be considered. 
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